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Abstract 

For nearly a decade, it lias been a mystery why the small average number of photons absorbed 
per atom from an ultrashort laser pulse is able to induce a strong magnetization within a few 
hundred femtoseconds. Here we resolve this mystery by directly computing the number of photons 
per atom layer by layer as the light wave propagates inside the sample. We find that for all the 24 
experiments considered here, each atom has more than one photon. The so-called photon shortage 
does not exist. By plotting the relative demagnetization change versus the number of photons 
absorbed per atom, we show that depending on the experimental condition, 0.1 photon can induce 
about 4% to 72% spin moment change. Our perturbation theory reveals that the demagnetization 
depends linearly on the amplitude of laser field. In addition, we find that the transition frequency 
of a sample may also play a role in magnetization processes. As far as the intensity is not zero, 
the intensity of the laser field only affects the matching range of the transition frequencies, but not 
whether the demagnetization can happen or not. 

PACS numbers: 75.40. Gb, 78.20.Ls, 75.70.-1, 78.47.J- 
Kcy words: 
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I. INTRODUCTION 



The pioneering discovery by Beaurepaire et al. [ll, that a femtosecond laser pulse can 
demagnetize Ni on a subpicosecond time scale or femtomagnetism, has inspired enormous 
scientific activities both experimentally and theoretically |2|-|7|. The significance of this 
discovery is that it demonstrates a possibility for nonthermal writing in a ferromagnetic 
medium. One prominent example is the inverse Faraday effect, where the laser can non- 
thermally switch spins js]. This process works even better when the temperature is lowered 
9|. Its potential applications require a good understanding of the underlying excitation 
mechanism. In spite of extensive investigations in this field 10, ll|, how the light transfers 
the photon energy to the system and subsequently demagnetizes the sample is still puzzling 
though new theoretical and experimental investigations emerge 3l- At the center of the 
debate is whether there are enough photons absorbed per atom (estimated at 0.01 12|) 
for magnetic moment change js], [12I2II. On the one hand, Koopmans et al. argue that 
the excitation density is too low to induce any substantial change in magnetization 13l-ll5|. 
If the number of photons is not enough in the first place, the excitation density must be 
very low. On the other hand, nearly all the experiments report the magnetization change 
either directly or indirectly. A strong demagnetization is incompatible with a shortage of 



photons. This puzzle affects our confidence in femtomagnetism |22l. |23|. Therefore, resolving 
this apparent contradiction is of paramount importance to femtomagnetism and its future 
applications. 



15 



2^ stated 



To this end, there are very few detailed investigations. Koopmans et al. 
that an effective photon number for demagnetization can only quench the magnetic moment 
of 10^"^ yUB/atom, which is much lower than the observed demagnetization of 0.003 yUe/atom. 
In 2007, Dalla et al. 1^ investigated the influence of photon angular momentum on the 
ultrafast demagnetization in Ni. Their results excluded direct transfer of angular momentum 
to be relevant for the demagnetization process and showed that the photon contribution to 
demagnetization is less than 0.01%. This motivated them and others to search for alternative 
mechanisms for strong demagnetization besides the spin-orbit coupling based mechanism 



25 1 . Such argument is reiterated by Stanciu et al. in Refs. [8|, ll2|]. Very recently, Hertel 



22l |. in a viewpoint on Ref. 



23|, claimed that the apparently simple assumption of a direct 



transfer of the photon spin to the magnetic system is not the solution. Majority of research 
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merely avoid this contradiction by referring to the theoretica 



work 



15| and the circularly polarized experiment in nickel 



argument made in the early 



previous paper 



16| . As we pointed out in a 



171] , the insensitivity of magnetization change to the light polarization is not 



a sufficient condition to rule out the direct involvement of photons in the first place. Laser 
affects the spin moment change in two ways. One is that the light changes the magnetic 
angular momentum m^. This is the case for the circularly polarized light. The other way is 
that the light changes the angular momentum /. This is the case for both the circularly and 
linearly polarized lights. If the magnetization change is not sensitive to circularly polarized 
lights, it only means that the channel is not effective, but it can not exclude the /- 
channel. As a result, one can not exclude the photon mechanism. A very recent theoretical 



investigation by Woodford 



181 ] reinforces this concept. Nevertheless it is important to note 



that irrespective of underlying mechanisms, such a low photon number is unlikely to induce 
a substantial magnetization in a sample. 

In this paper, we develop a generic scheme to compute the average photons absorbed 
per atom and show that for all 24 sets of experimental data considered here, each atom has 
more than 1 photon. For a weak laser field, we examine the relation between the relative 
demagnetization change and the mean photons absorbed per atom. Our results show that the 
small number of photons absorbed per atom can induce a strong magnetization. Moreover, 
the linear dependence of demagnetization change on the amplitude of laser field is consistent 
with our perturbation result. For the strong laser intensity, we resort to the two-level model 
system. We find that an effective demagnetization change occurs even with a weak laser 
field as far as the system is at resonance. 

This paper is structured as follows. Section II presents a formal algorithm to compute 
the photon number in femtosecond magnetism. In Sec. Ill we show the relation between the 
relative demagnetization change and the mean photons absorbed per atom. Then we present 
a theoretical investigation on the demagnetization change from the perturbation theory in 
Sec. IV. For a strong field, results are presented in a two-level model in Sec. V. Finally, the 
main conclusions of our study are summarized in Sec. VI 
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II. PHOTON NUMBER IN FEMTOSECOND MAGNETISM 



We start with a laser pulse propagating along the positive z axis with the field along the 

% clXlS, 

£^a;(c<;, r, 2;; t) = exp (— t^/r^) exp iuj{^ — t^ (1) 

where is the amplitude of laser field, r is the pulse duration, v is the phase velocity, uj 
is the laser frequency and t is the time. Since v = c/{n + ik), we can rewrite the above 
equation as 



/nz 



Sx{uj,T, z;t) = Aoexp {—t /t ) exp iuj i t) exp ( z) (2) 



V c 



ujk 



c 



where n and k are real and imaginary parts of the index of refraction, respectively, and both 
are wavelength-dependent, c is the speed of light. 

The laser intensity / for the linearly polarized light is computed from 26 1 

I{uj, r, z; t) = 2neoc\Sx{io, t; t)p = 2neoc/lQexp (— 2t^/r^) exp ( z) 

2; 

= /o(u;,r;t)exp(--) (3) 

where eo is the permittivity of free space and /o(w, t; t) is the laser intensity before penetrat- 
ing the sample. This is the well-known Beer-Lambert law. Here d is the penetration depth, 
defined a.s d = c/2uk = X/Airk, which describes how the light intensity falls off starting 
from the surface of a sample. It is clear that d itself is intensity- independent. Since we are 
interested in the pulse energy fluence F{u, r, z), we integrate the above equation over time 
and find, 

/oo 
r, 2;; (4) 
-oo 

which can be simplified as 

/TT Z Z 

-rexp (--) = Fma^exp (--). (5) 

where F^ax is the initial laser fiuence given in experiments. This is the exact expression for 
the laser energy fiuence at depth z, if the pulse is a Gaussian function. 

If a laser pulse of fiuence F{u, r, z) shines on a spot with area A, the total number of 
photons within A at depth z is 

Nphoton{(^, z) = -^^j— exp (-^) (6) 
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where hu is the photon energy, and fi is the Planck's constant over 27r. This equation re- 
veals some crucial information: (i) Nphoton is a surface quantity; (ii) the number of atoms 
illuminated by these photons must be a surface quantity as well. In other words, we must 
compute the photon numbers provided per atom layer by layer, since the light wave propa- 
gates uniaxially. Consider an fee structure with lattice constant a and surface area A, the 
number of atoms in each layer is 

2A 

Natom = — ^ (7) 

where 2 comes from the fact that each unit cell has two atoms per layer. The mean number 
of photons available to each atom at different depth z is 

where S^ax is the maximum mean number of photons provided to each atom in the top 
layer. This equation gives the true mean number of photons provided to each atom for 
demagnetization at depth z. 

Under the low pulse energy fiuence, ignoring small reflected photons, which is justified in 
highly absorbed metals, we estimate the photons absorbed per atom at the jth layer as 



J = l,2,3,--- (9) 



To appreciate how large S and r/j are, in Table I we list the results from 24 different 
sets of experimental data increasing from top to bottom. This table is very telling. The 
maximum photon number 'B.^ax ranges from 1.368 to 90.430. The photon number at the 
penetration depth varies from 0.503 to 33.267. Note that since majority of samples are 
thinner than the penetration depth, the number of photons available to each atom exceeds 
1. Here we emphasize that there are enough photons for each atom. Absorbed photon in the 
first layer 771 ranges from 0.017 to 1.075. The mean photons absorbed per atom within the 
sample is less than r\i. Koopmans' group [l^ has the smallest value of 0.017. Beaurepaire 
et al. [l|'s rji value is 0.176, ten times higher than that of Koopmans' rji. Cheskis et al. 
{35! have the largest rji value of 1.075 or roughly one photon per atom, and importantly 
they already found that the demagnetization is saturated. This provides a first indication 
that one does not need a large amount of photons to demagnetize the sample. The small 
number of photons absorbed does not necessarily mean that they can not induce a strong 
magnetization change. 
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III. DEMAGNETIZATION CHANGE VS PHOTONS ABSORBED 



We need to map out the connection between the magnetization change and number of 
the photons absorbed. We investigate the relative demagnetization changes AM/M which 



are extracted from their corresponding 



a weak laser intensity 



[il,0,[3,Q[3: 



iteratures versus the photons absorbed per atom for 



20 



35| . We purposely choose a weaker intensity in 



order to search for a one-to-one correspondence between the number of photons absorbed 
and the amount of magnetization change. For a strong laser, due to the saturation, such a 
relation will become complicated. To find some valuable information on the demagnetization 
mechanism associated with the number of photons, it's necessary to estimate the mean 
photons absorbed per atom f] more accurately. Based on Eq. (9), f] is computed from 

1 

^ = ^E^^- (10) 

where k is the number of layers within the sample's thickness. If the samples' thickness is 
thinner than the penetration depth, we calculate f] averaged within the sample's thickness, or 
else within the penetration depth. The f] itself includes the reflected and absorbed photons. 

Figure 1 shows AM/M as a function of f/. Four solid lines (guide for the eye) from top 
;o bottom represent the results by Bigot et al. Beaurepaire et al. jlj], Cheskis et al. 



35 
13 



and Atxitia et al. (20 
and Wilks et al 



, respectively. The corresponding results from Koopmans et al. 



19| are also shown in the bottom left corner of the figure. The 
vertical dotted line shows that 0.1 photon can induce about 4% to 72% spin momentum 
change for different laser durations. This figure is very insightful, (i) The slope of these 
lines approximately characterizes the demagnetization ability of a laser pulse. Each solid line 
denotes the results obtained by an identical laser pulse which has the same pulse duration; 
(ii) The amplitude of laser Aq plays a key role in the demagnetization change. For the 
same fluence, as predicted by Eq. (5), the smaller the pulse duration is, the bigger the 
laser amplitude becomes. Within the dipole approximation, the interaction between the 
laser and the system is dominated by the amplitude of laser. Therefore, for same f/, the 
amplitudes increase from bottom to top; (iii) With the same pulse duration parameter, the 
demagnetization change increases along with the photons absorbed per atom or the laser 
fluence. The solid lines with slopes of 1.5, 3.8 and 7.4 indicate that one does not need a large 



amount of photons to induce a substantial moment change. For Koopmans et al. 



13 
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and Wilks et al. 



19| data, because their laser energy fluences are very small and the pulse 



durations are long, the demagnetization changes are much smaller, but are still consistent 
with the above picture. 

From Fig. 1, it is obvious that the laser amplitude plays a crucial role in the demagneti- 
zation changes. However, is the laser amplitude the only deciding factor for demagnetization 
change? The answer is negative. In fact, using the concept of photons absorbed per atom 
to explain the strong demagnetization change obviously neglects two very important fac- 
tors. First, it doesn't take into account the interaction between the light and the material. 
Light has the energy (how strongly the field oscillates) and the frequency (how fast the field 
oscillates with time), fj only takes into account the energy but not the frequency, nor the 
transition matrix elements between different states. Second, once the number of photons 
per atom becomes small, it is well known that the quasi-classical description of photons 
absorbed per atom for demagnetization becomes invalid. In particular, when the average 
number of photons provided to each atom is less than 40, the electric field behaves quantum 
mechanically 38|, i.e., the field oscillates strongly around its average value. The less the 
photons are per atom, the stronger the oscillation is. Therefore, it is necessary to reconsider 
this hurdle from a different perspective, and we examine this issue from the perturbation 
theory and two-level model, respectively. 



IV. PERTURBATION THEORY AND WEAK INTENSITY 

In Sec. Ill, it reveals that the relative demagnetization change AM/M is proportional to 
the mean photons absorbed per atom f/ under a weak laser field. This allows us to treat the 
laser field perturbatively. The Hamiltonian of a system can be described hy H = + , 
where is the time-independent Hamiltonian of the unperturbed system, = —ejl- £{t) 
is the time-dependent perturbation. We start with the Liouville equation ihp = [H, p] for 
the density matrix. We keep only the first order term and have 

thp^'^ = [H',p^'^] + [H',p^% (11) 

where p^^^ and p^^^ are the zeroth and first order of the density matrices, respectively. If we 
make an unitary transformation as p^ 

1) = ^-iH^t/nQ^iH^t/h^ g^^^ Eq. (11) can be written as 
^/^g = e^^^°*/^[i^^pW]e-*^°*/^ (12) 
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We integrate this equation over time and take oo) = 0. Then apply the eigenstate {n\ 

on the left and \m) on the right of this equation, and we can obtain 

P^lit) = Up^:1 - P^±)ePnme-'--'^^' f dt'£{ty--''-^'' (13) 

where Unm = {En — -E'm)/^? Pnm the transition matrix elements, and F the damping factor. 
Theoretically, the demagnetization change is involved in the time-dependent density matrix 
p'nlit) through MP{t) = Tr[S,p^^\t)], where S, is the spin matrix. By investigating the 
density matrix, we can reveal some crucial details of the demagnetization. Next we discuss 
two typical cases. 

Case 1, Continuum wave laser. Consider a periodical field perturbation = — e/xAoe"*'^*, 
the density matrix becomes 

1 _ e/inm, (0) (0) ^ A _2rt ( siu - T COS wt) + ( Aw COS + T siu a;t)z 

PnmV') — KPnn ~ PmrnJ^OG Acj^ + 

where Au = u — Unm- Note that its frequency dependence is exactly same even if we treat 
photons quantum mechanically. It shows that the density matrix pnm{t) not only depends on 
the amplitude of laser field Aq which is consistent with the conclusion obtained from Fig. 1, 
but also depends on the resonance term ( Aa;^Vr^ ) • Since the spin moment change is directly 
proportional to pnrn{t), this predicts that the magnetic moment change will depend linearly 
on the amplitude of laser field in the weak field limit. On the other hand, the resonance term 
demonstrates that even for the exactly same laser amplitude, the transition or the magnetic 
change can be very different for different frequencies. This has been largely ignored in the 
literature. 

Case 2, Pulse laser. We assume that laser field has a form = — e/iy4oe~*^'''^^e~*'^*. Eq. 
(13) becomes 

P^:iit) = '-^{p^^l-p^±)e—'~''Ao f rft'e-'^/^^e-^^-'-^*' 

^'-^{p^l-p^±)e—'-''R{t) 
in 

where R{t) is defined as the module of response function. To get some insightful information 
about the perturbation of laser field, we choose r = 12 fs, and Aq = 0.05 V/A, and integrate 
the response function numerically. 

Figure 2(a) shows the module of response function R{t) as a function of time for twelve 
different energy detuning AE {TiAoj) from 0.05 to 0.21 eV. We can see that as time goes by, 
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R{t) first increases to a peak value, and then settles down to its final value. Importantly 
the peak values of R{t) vary a lot for different AE even with the same laser field amplitude. 
The biggest one is about 0.9 {AE = 0.05 eV) and the smallest one is about 0.07 {AE = 
0.15 eV). Their difference is over an order of magnitude. It demonstrates that the transition 
probability or spin momentum change can be very large when AE becomes very small or 
the system is excited resonantly, even if the field intensity is weak. 

Figure 2(b) compares the peak and final values of R{t) as a function of AE. The final 
R{t) values decrease more quickly with AE than that of peak ones. The maximum difference 
is 0.16 at AE = 0.22 eV, and the minimum difference is 0.09 at AE = 0.4 eV. When AE 
reaches 0.07 eV, the peak and final values are almost same. It implies that regardless of 
the amplitude of laser field the transition from one state to another state is finished when 
AE goes to zero. Figure 2(c) shows the time needed reaching the peak value as a function 
of AE. This indicates that the peak time is reduced as AE becomes larger. All the above 
results are obtained within a weak field limit. Once the laser field becomes strong, the first 
order perturbation becomes invalid. Under this situation, we resort to a two-level model 
system. 



V. TWO-LEVEL MODEL AND STRONG INTENSITY 



The two-level model has been extensively used in atoms 39|, |40|], semiconductors 41 



and ferromagnetic materials 4J]. For a system excited by a laser, the two-level model can 
give us a quantitative understanding of the spin transition. Here we directly quote the Eq. 
(7) of Ref. 4^ in our previous work about the spin change ASz for a transition from state 

I a) to state \b) 

AS, = Kb{t;oo){{b\s,\b)-{a\sM) (16) 

where t is time, and u is laser frequency. ^ab{t] ou), the probability amplitude of finding the 
system at time t in state \b), is equal to 



Kb{t;uj) 





Wab 


2 







sm 



Wab? 



+ {bJ - Uba) 



where Wab = fJ^ab^o- The transition matrix elements fiab can be obtained 



corresponding momentum operator Pab from the ab initio calculation 



(17) 



3y calculating the 



45N47| 
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The probability u) is an oscillatory function of time; for certain values of t, 

ii.abit] uj) = 0, meaning that the system returns to the initial state \a). At the same time it 
also reveals fast oscillation expressed by the latter term of Eq. (17). If \Wab\ is large enough, 
the value of Ka;,(t; u) or AS^ can be very large. For the off-resonance excitation, the laser 
field amplitude Aq controls the demagnetization change entirely. This happens in semicon- 
ductors. For the resonance excitation, or u = Uba, regardless of how weak the perturbation 
is, the field can cause the system to transit from state \a) to state \b) 48|. Importantly, 
the field intensity only affects the time needed for the system to transit from \a) to \b), not 
whether a transition can occur or not. The smaller the field intensity is, the longer the time 
it takes. 

For a ferromagnetic metal like nickel, the chance that the transition frequency u^a matches 
that of the laser field is very high 49|] . This explains the observed demagnetization in those 



experiments, in spite of a relatively weak laser electric field. Figure 3 shows the detailed 
dependence of the spin change on the laser electric field. For a weak laser field intensity 
with 0.001 V/A, the spin change is possible as far as Uba matches u. This demonstrates 
that the strong demagnetization change is achievable in experiments even if the laser field 
intensity is very weak. If the laser field becomes larger, the range of matching frequency 
becomes broad (see Fig. 3). In semiconductors, Pavlov et al. 50|] showed that increasing 
temperature reduces the GaAs band gap. This may be a test case for our theory. 



VI. CONCLUSION 



We have clarified a long-standing conceptual puzzle of the photon shortage in femtosecond 
magnetism by comparing the relative demagnetization change versus the mean photons 
absorbed per atom. In the weak laser field, it increases along with the mean photons 
absorbed per atom f/, which is consistent with the perturbation theory. Importantly, the 
results show that few photons absorbed per atom can induce considerable spin moment 
changes. The spin moment can be reduced at resonance even if the field is weak. Our 
findings overcome a big hurdle in ferromagnetism and should inspire new experimental and 
theoretical investigations into the role of photons and their interactions with electrons in 
the magnetization change. 
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TABLE I: The calculated number of photons per atom in Ni. H^^^. is the maximum photon number 
per atom, d is the optical penetration depth. is the photon number at d. r]i is the photons 

absorbed per atom in the first top layer. The other parameters are taken from the literature, n 
and k are real and imaginary parts of the index of refraction , respectively. 



No. 


F 

(mJ/cm^) 


A 

(nm) 


n + ik 


d 
(nm) 






'/I 


Ref. 


1 


0.6 


729 


2.28+4.18f 


13.93 


1.368 


0.503 


0.017 


[141 


2 


0.76 


729 


2.28+4.18f 


13.93 


1.741 


0.641 


0.022 


[131 


3 


1.4 


790 


2.46+4.35i 


14.45 


3.457 


1.272 


0.042 


[191 


4 


1.45 


798 


2.476+4.375i 14.52 


3.617 


1.331 


0.044 


[71 


5 


1.8 


800 


2.48+4.38^ 


1 A A r 

14.45 


4.501 


1.656 


0.055 


[281 


6 


2.04 


729 


2.28+4. 18^ 


13.93 


4.644 


1.708 


0.058 


[131 


7 


2.0 


785 


2.45+4.34? 


14.50 


4.908 


1.805 


0.059 


[161 


8 


2.0 


800 


2.48+4.38^ 


14.45 


5.001 


1.840 


0.061 


[291 


9 


2.5 


800 


r\ A a \ A O O ' 

2.48+4.38? 


1 A A r 

14.45 


6.252 


2.300 


0.076 


[Hi 


10 


2.83 


800 


r\ A ci 1 A O O ■ 

2.48+4.38^ 


1 A A r 

14.45 


7.075 


2.603 


0.086 


[301 


11 


3.5 


800 


9 48-1-4 


1 4 d'^ 


8.752 


3.220 


0.106 


111] 


12 


4.4 


800 


2.48+4.38i 


14.45 


11.003 


4.048 


0.133 


[Hi 


13 


5.3 


800 


2.48+4.38i 


14.45 


13.254 


4.876 


0.161 


[Hi 


14 


7.0 


620 


1.93+3.65i 


13.53 


13.566 


4.991 


0.176 


[1 311 


15 


6.0 


800 


2.48+4.38i 


14.45 


15.004 


5.520 


0.182 


[321 


16 


8.0 


620 


1.93+3.65i 


13.52 


15.504 


5.704 


0.201 


[311 


17 


7.1 


800 


2.48+4.38i 


14.45 


17.755 


6.532 


0.215 


[28] 


18 


12.0 


800 


2.48+4.38i 


14.45 


30.008 11.039 0.364 


[33] 


19 


12.73 


827 


2.53+4.47f 


14.73 


32.897 12.102 0.391 


[10] 


20 


13 


827 


2.53+4.47i 


14.73 


33.588 12.356 0.399 


m 


21 


13.3 


827 


2.48+4.38i 


14.73 


34.363 12.642 0.409 


[35] 


22 


15.72 


800 


2.48+4.38i 


14.45 


39.308 14.461 0.476 


[36] 


23 


20 


800 


2.48+4.38i 


14.45 


50.014 18.400 0.606 


[33 


24 


35 


827 


2.13+4.17i 


14.73 


90.430 33.267 1.075 


[35] 
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ri (Photons/atom) 

FIG. 1: Dependence of the relative demagnetization changes on the average photons absorbed per 
atom for six different sets of experimental data. The solid lines are guides for the eye, whose slopes 
are shown on the lines. The laser pulse durations and the photon wavelengths are also given near 
their data. 
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FIG. 2: (a) The module of response function R(t) as a function of time at different AE. (b) The 
dependence of the peak and final values on AE. (c) The time needed to reach the peak as a 
function of AE. 
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FIG. 3: Amplitude of spin momentum change as a function of the detuning AE = fi{oj — ujha) 
for three groups at different laser electric fields. The amplitude of laser field is deduced from the 
experimental laser fluence based on the Eq. (5). All the curves are shifted to zero at resonance. 
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